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Abstract 

We show how one can measure anomalous WW Z- and WW^- 
couphngs with minimal statistical error using integrated observables, 
without having to assume that the anomalous couplings are small. 
We propose a parametrisation of these couplings which is well suited 
for the extraction of both single and many parameters, and which 
leads to a very simple form of the integrated cross section, from which 
additional information on the couplings can be obtained. 
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1 Introduction 

The direct and precise measurement of the self-couphng between the elec- 
troweak gauge bosons in VT-pair production will be a crucial step in testing 
the standard model of electroweak interactions and searching for physics be- 
yond it. It will form an important part of the physics programme at LEP2 
and at a planned linear e"'"e~-collider (LC). As is well known there are three 
diagrams at tree level that contribute to the amplitude of e^e~ — ^ W^W~ 
in the standard model, one with t-channel neutrino exchange and the other 
two with a 7 or Z in the s-channel, involving the vertices WW'-f and WW Z. 
One can parametrise the corresponding vertex functions in order to quantify 
the couplings and to compare them with their form in the standard model. 
In the most general form respecting Lorentz covariance each vertex involves 
seven complex form factors [|^], three of which give couplings that violate CP 
symmetry. 

Without further physical assumptions one is thus left with 28 real param- 
eters whose simultaneous extraction in one experiment looks quite hopeless. 
Given the limited event statistics expected at both LEP2 and the LC one will 
only obtain meaningful errors on a reduced number of coupling parameters 
at one time. This may be achieved by imposing certain constraints on the 
full set of coupling constants; various suggestions for such constraints based 
on symmetry considerations have been made in the literature P, §|. One 
must however keep in mind that experimental values or bounds on couplings 
that have been obtained with particular constraints cannot be converted into 
results without constraints or with different ones; the information lost by as- 
suming relations between couplings cannot be retrieved. Although imposing 
such constraints is certainly legitimate and can be useful we stress that a 
data analysis with independent couplings will be valuable, both from the 
point of view of model independence and the capability to compare results 
of different experiments. 

We remark that of course one can also give (reasonably small) errors on 
single or few couplings in a multi-parameter analysis. In this paper we pro- 
pose a parametrisation of the couplings which is well adapted to this end, the 
statistical errors on the different measured parameters being approximately 
uncorrelated. 

We will work in the framework of optimal observables, a way to extract 
unknown coupling parameters introduced for the case of one parameter in 
IQ, that has since been used for various reactions [^, 0; ^- General aspects of 
this method, in particular its extension to an arbitrary number of parameters, 
as well as its application to W^W~ production were discussed in |Q. In this 
paper we investigate again the reaction e'^e'^ -^ W^W~ . We concentrate 



here on the decay channels, where one W decays hadronically and the other 
into an electron or muon and its neutrino. Calculated with the Born level 
cross section of the standard model the statistics of these channels is about 
3000 events for a collision energy of ^/s = 190 GeV and 500 pb~^ integrated 
luminosity, which are typical planned LEP2 parameters, and about 22000 
events with 10 fb~ at ^/s = 500 GeV, which might be achieved at the LC 
A complementary source of information is the integrated cross section, 
which is a quadratic function of the triple gauge couplings. The combination 
of information from the total event rate and from observables that make use 
of the detailed distribution in the final state has for example been used in 



T0| , where CP violation in the decay Z -^ bbg was investigated. 

In sec. ^ of this paper we will further develop some aspects of the method 
of optimal observables, in particular we will show how to apply it without 
the linear approximation in the coupling parameters that was used in 0. In 
sec. 0we then propose a parametrisation of the couplings that simultaneously 
diagonalises certain matrices connected with our observables and with the 
integrated cross section. These parameters achieve two goals: their quadratic 
contribution to the total cross section is a simple sum of squares and the co- 
variance matrix of the corresponding optimal observables is diagonal. The 
methods which we use for this purpose are borrowed from the theory of small 



oscillations of a system with / degrees of freedom (cf. e.g. |TT|). Our param- 
eters correspond to "normal coordinates" and their use in an experimental 
analysis should in our view present several advantages. We give some numer- 
ical examples for ly-pair production at LEP2 and the LC in sec. ^ and make 
some further remarks on how our proposal might be implemented in practice 
in sec. ^. The last section of this paper gives a summary of our main points. 

2 Optimal observables: analysis beyond lead- 
ing order 

The method of optimal observables has previously been presented in the 
approximation that the couplings to be extracted are sufficiently small to 
allow for a leading order Taylor expansion of various expressions. Here we 
show how to use it beyond this approximation. 

Let us denote by gi the real and imaginary parts of the WW'j and WWZ 
form factors minus their values in the standard model at tree level. As 
the amplitude of our process is linear in these couplings we can write the 



differential cross section as 

— = So{(j)) + J2 SiA<t^) 9i + Yl '52,ij(</>) 9i9j , (1) 

wfiere S2,ij{4>) is a positive semidefinite symmetric matrix, (j) collectively 
denotes the set of measured phase space variables. The integrated cross 
section is 

a = (Toil + Y, ^i,i 9i + Y. '^2,2i 9i9j (2) 

\ i ij J 

with the standard model cross section ctq = J d(j) Sq{(J)) and coefficients 

d-i,i = — d(j)SiA(j)) ' cr2,ij = — d(j) S2,ij{(j)) . (3) 

aoJ (JqJ 

The idea of using integrated observables is to define suitable functions Ci(0) 
of the phase space variables and to extract the unknown couplings from their 
measured mean values Oi. Let us give the details. From (|l]) and @) we 
obtain the expectation value E[Oi\ of Oi as 

H Cij 9j + J2 (lijk 9j9k 

E[0,] - E,[0,] = --^^ f— .^ (4) 

1 + Z^ ^1 ,i 9j + 2^ ^2,jfc 9j9k 

j jk 

with the standard model expectation value -E'olC'j] = {J d(f)OiSo)/(TQ and 
coefficients 



Qijk 



- fd(j)OA,,-Eo[0,]a^,, , 

(To J 

- fd(l)OiS2,,k-Eo[0,]a2,jk ■ (5) 

CTn J 



We remark in passing that the coefficients in (^ can be written in a compact 
form as 

Cij = VolOi , Sij/So] , Qijk = Vo[Oi , S2jk/So] , 

o"ij = Eo[Sij/So] , (T2,jfc = Eo[S2jk/So] , (6) 

where Vo[f, g] = Eo[fg] —Eo[f] Eo[g] is the covariance of /(0) and g{(j)) in the 
standard model. Note that (T2,jfc is symmetric and positive definite, whereas 
Qijk as a matrix in j and k is symmetric but in general indefinite. 



An estimation of the couplings can now be obtained by solving the system 
@ with E[Oi] replaced by the mean values Oi, 

Y. Cii 9j + Y. %k 9j9k 

O. - Eo[0.] = --^^ j^^ , (7) 

1 + 2^ ^i,i 9j + 2^ ^2,jfc 9j9k 

provided of course one has n observables for n unknown couplings. When 
the system (|^) is linearised in the gi it is easily solved by inversion of the 
matrix Cij. One is however not constrained to do so and can instead solve 
the exact set of equations (|^. By multiplication with the denominator it 
can be rearranged to a coupled set of quadratic equations in the gi and will 
in general have several solutions. Some of these may be complex and thus 
ruled out, but from the information of the Oi alone one cannot tell which of 
the remaining real ones is the physical solution. We will come back to this 
point. 

The measured mean values Oi are of course only equal to the E[Oi] up to 
systematic and statistical errors. We only consider the latter here, which are 
given by the covariance matrix V{0)ij of the observables Oi divided by the 
number A^ of events in the analysis. To convert the errors on the observables 
into errors on the extracted couplings we use the quantity 

x\g)=Y{(^^-E[0,])NV{Or\,{0,-E[0,]) , (8) 

which depends on the gi through the E[Oi] given in (Q). Solving (|^) is 
tantamount to minimising x^ with Xmin — O5 ^^"^ ^ confidence region on the 
couplings is as usual given by 

X^{9) - xliin < const. (9) 

with the constant determined by the desired confidence level. 

There are several possible choices for the covariance matrix V{0) in (|]). 
It can be 

1. determined from the measured distribution of the observables Oi, 

2. calculated from the differential cross section (0), taking for the gi the 
values extracted in the measurement, 

3. calculated for vanishing couplings gi, 

4. calculated as a function of the couplings. 



Choices 1. and 2. should lead to the same results in the limit of large A^ 
where the statistical errors on the measured V{0)ij and Qi become small. 
Comparison of the covariance matrices obtained by these two methods might 
indeed be helpful to rule out unphysical solutions of (|^). Choice 3. in turn 
will be a good approximation of 2. if the couplings are small enough. We 
consider possibility 4. as the least practical one, except maybe for the case of 
one coupling. For several couplings the expression of V{0)ij as a function of 
the Qi involves tensors of rank up to four and is even more complicated than 
the one for the expectation values (^, and the inverse matrix is yet more 
clumsy. For this reason we will discard choice 4. in the following. 

In we considered an analysis at leading order in the gi, where one uses 
the linearised form of (^ to estimate the couplings: 

gj = T.^~'jk{c>k-Eo[Ok]) . (10) 

k 

Correspondingly the linear approximation of (^ is used in the expression (H) 
of x^ which then reads 

x\9) = T.i9^ - 9^) V{g)-\, [g, - g,) , (11) 

where 

V{g)-^ = Nc^ -V^O)-^ -c (12) 

is the inverse covariance matrix of the estimated couplings [^ . As one works 
to leading order in the gi one can approximate V{0) by its value for zero 
couplings, i.e. choose possibility 3. above. 

The confidence regions x^ig) < const, for the measured couplings are 
then ellipsoids in the space of the gi with centre at {gi). The optimal observ- 
ables 

O.W = ^-^ (13) 

discussed in [^ have the property that to leading order the statistical errors 
on the estimated couplings are the smallest possible ones that can be ob- 
tained with any method, including e.g. a maximum likelihood fit to the full 
distribution of given by the differential cross section (p. 

Note that one can still use the linearised expressions ([ToD and (pUD in 
an analysis beyond leading order. The error x^id) < 1 on the couplings 
will be given by an ellipsoid with defining matrix ([T^) , where V{0) is the 
covariance matrix at the actual values of the couplings. These errors will in 
general no longer be optimal, so that when the leading order approximation is 
not good one might obtain better errors with a different choice of observables. 



More importantly, however, the extracted values of the couplings are biased: 
averaged over a large number of experiments the measured couplings differ 
from the actual ones by terms quadratic in the gi. 

If instead one uses the full expressions (H), (^ and (H) one has no bias 
on the extracted coupling parameters, provided the number A^ of events in 
the analysis is large enough. Let us see if we can find optimal observables for 
this case. To this end we expand the differential cross section around some 
values (ji of the couphngs: 

-T- = So + J2 Si,i{<P) ■ i9i -9i) + J2 ^2,ij{4>) ■ {9i - 9i){9j - 9j) ■ (14) 
f i ij 

The corresponding zeroth order cross sections and mean values are ao = 
J d(f) So and -Eo[C^j] = (/'^0C'j5'o)/cfo, respectively. We then can re-express 
E[Oi] in (^, replacing gi with g^ — gi, Eq with Eq, and using new coefficients 
Cij etc. constructed as in (^, (|^). Making the same replacements in (J^) we 
have an alternative set of equations to extract the coupling parameters. 

It can be shown that for sufficiently large N the confidence regions ob- 
tained from d^), (^ in a nonlinear analysis are again ellipsoids given by 
X^{9) ^ const. One can then write x^{9) ^^ iii (0)) but with V{g)^^ of ([T2|) 
replaced by 

Vig)-^ = Nd'' -ViOy^-c , (15) 

where c corresponds to an expansion (0) of da/d(j) about the actual values 
of the couplings. The main point of the argument is that for large N the 
statistical errors on the Oi become small, so that the extracted couplings will 
be sufficiently close to the actual ones to allow for a linearisation of (^) and 
(0), cf. ig, p.695, and [|T3]. 



Finally one can construct new observables 

om = ^ (16) 

from (|14[). They will be optimal, i.e. have minimum statistical error if the (ji 
are equal to the actual values of the g^. In the appendix we show that, up 



to linear reparametrisations given in (39), this is the only set of n integrated 



observables that measures the n couplings with minimum error. There is 
hence no choice of observables that would be optimal for all values of the 
actual coupling parameters. As these are unknown one can in practice not 
write down the truly "optimal" observables, but our argument tells us how 
one can improve on the choice in (0) if one has some previous estimates gi of 
the couplings (cf. also 0). One may then choose to perform a leading order 
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analysis as described above, linearising about gi = gi instead of gi = Q. A 
practical way to proceed could be to estimate the parameters gi at first using 
the linearised method around gi = 0. Suppose this gives as best estimate 
some values c/i. Then in a second step one could set gi = g^ and use the 
linearised method around gi to improve the estimate etc. 

At this point we wish to comment on the "optimal technique" for deter- 
mining unknown parameters in the differential cross section that has been 
proposed in [Q. The "weighting functions" Wi{(j)) there depend on the actual 



values of the parameters one wants to extract and are thus not "observables" . 
Only if one sets the unknown parameters in the Wi{(f)) equal to some previous 
estimates of them can one use these functions to weight individual events; 
the better these estimates are the more sensitive the functions will be. If one 



does this then the set Wi{(f)) is equivalent to our observables (|T6D defined for 
some estimates gi of the coupling parameters. 

We finally remark that if N is not large enough the statistical errors on 
the mean values Oi and thus on the measured couplings might be so large 
that they lead into a region where a linearisation of (Q) is not a good ap- 
proximation. The covariance matrix V{g) is then no longer given by (plsl). 
Moreover the errors on the couplings might be asymmetric and the shape 
of the confidence region defined by (|]), (|^) very different from an ellipsoid, 
so that knowledge of V{g) is not sufficient to estimate the errors on the gi. 
In such a case we cannot say on general grounds how sensitive our observ- 
ables are. Incidentally this also holds for other extraction methods such as 
maximum likelihood fits, whose optimal properties are realised in the limit 
N ^ oo. If one is rather far from this limit the sensitivity of a method 
will have to be determined by other means, e.g. by detailed Monte Carlo 
simulations. 

The method we have outlined can of course also be applied if one chooses 
to reduce the number of unknown parameters by imposing certain linear 
constraints on the couplings. One may still use the observables ( [T^ ) corre- 
sponding to the full set of couplings but minimise x^ iii ® for the reduced 
set; in this case one can of course not take choice 2. for V{0). In general 
Xmin is then different from zero and its value indicates to which extent the 
particular constraints on the couplings are compatible with the data. If N 
is large enough Xmm follows in fact a x^-distribution with n — m degrees of 
freedom for n observables and m independent couplings so that its value can 
be converted into a confidence level |I2| . 

We conclude with a remark on the use of optimal observables in practice. 
A realistic data analysis will not be good enough if the Born approximation 
of the differential cross section (||) is used. Both higher-order theoretical 
corrections, such as initial state radiation and the finite W width, and ex- 



perimental effects like detection efficiency and resolution will modify the 
observed distribution of the phase space parameters cf). If they are taken 
into account in the determination of the coefficients in (^, (|^ and of the 
covariance matrix V{0) they will not lead to any bias in the extraction of 
the couplings and their errors. While this will presumably be done with sets 
of generated events and might be computationally intensive one still has to 
determine only a rather limited number of "sensitivity" constants. On the 
other hand one needs to know the observables Cj(0) of ( ]T3[ ) as functions over 
the entire experimental phase space, so that the expressions of 5*0 and Si^i 
used to construct them will in practice be taken from a less sophisticated 
approximation to the actual distribution of (j) in order to keep them manage- 
able. The observables are then no longer optimal, and it will depend on the 
individual case which approximations of 5*0, Si,i are good enough to obtain 
observables with a sensitivity close to the optimal one. 

2.1 Discrete symmetries 

In it was shown how with a suitable combination of all semileptonic 
W~W^ decay channels one can define observables that are either even or odd 
under the discrete transformations CP and CPT, where C denotes charge 
conjugation, P the parity transformation, and T the "naive" time reversal 
operation which flips particle momenta and spins but does not interchange 
initial and final state. Under the conditions on the experimental setup and 
event selection spelt out in ||^ we have two important symmetry properties: 

1 . A CP odd observable can only have a nonzero expectation value if CP 
symmetry is violated in the reaction. 

2. If the expectation value of a CPT odd observable is nonzero the transi- 
tion amplitude must have an absorptive part whose phase must satisfy 
certain requirements in order to give an interference with the nonab- 
sorptive part of the amplitude. 

We assume in this analysis that any nonstandard physics in the reaction 
is due to the triple gauge vertices. In the standard model one needs at 
least two loops to violate CP; to a good accuracy the triple gauge couplings 
are therefore the only possible source of CP violation. For our process, 
i.e. e^e~ annihilation into four fermions, an absorptive part that satisfies 
the requirements mentioned in point 2. will appear in the Standard Model 
already at next-to-leading order in the electroweak fine structure constant, 
either through nonresonant diagrams or through loop corrections. To leading 



order, however, they are only due to the imaginary parts of triple gauge 
couplings. 

In this approximation the optimal observables (|T^ are CP even (odd) if 
they correspond to CP conserving (violating) couplings, and CPT even (odd) 
if they correspond to the real (imaginary) parts of form factors. The coeffi- 
cient matrix cij is then block diagonal in four symmetry classes of observables 
and three-boson-couplings: 

(a) : CP and CPf even 

(6) : CP even and CPf odd 

(c) : CP odd and CPf even 

{d) : CP and CPf odd. 

In the leading order analysis one thus can treat these four classes of cou- 
plings separately and benefit from a great reduction of unknown parameters. 
Beyond leading order, however, form factors of any symmetry can contribute 

to E[0;\: 

• in the integrated cross section and thus in the denominator of -E'[Oj] — 
Eq [Oi\ in (|^) couplings of all four classes enter quadratically, couplings 
of class (a) also appear linearly; 

• if Oi belongs to class (a) the numerator of E[Oi\ — -EolC'i] has terms 
linear in the couplings of this class but couplings of all four classes enter 
quadratically through qijk ; 

• if Oi belongs to a CP ( CPT) odd coupling then the numerator in (^ is 
only linear in CP ( CPT) odd couplings, but it contains also quadratic 
terms where a CP ( CPT) odd coupling is multiplied with a CP ( CPT) 
even one. 

We remark that this leads to different behaviours of -E[Oj] as one or more 
couplings Qi become large: whereas for observables in classes (6), (c) and (rf) 
the expectation value goes to zero when a coupling of the same class goes to 
plus or minus infinity the corresponding limit of an observable in class (a) 
can be a positive or negative constant or zero. 

In a nonlinear analysis one will therefore in principle have to consider 
couplings with all symmetries at the same time. In practice one might choose 
simpler procedures if the linear approximation is expected to be not too bad 
and if one wants to calculate corrections to it. One might for instance first 
analyse the four symmetry classes separately, neglecting in each case the 
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contributions of the three other classes at the r.h.s. of (|^ and then refine 
the analysis of a class by taking the values obtained in the first step for the 
couplings in the other classes as fixed in (|^. 

We emphasise that even beyond the leading order approximation it is 
still true that a nonzero mean value of a CP or a CPT odd observable is an 
unambiguous sign of CP violation or the presence of absorptive parts in the 
process, respectively. The extraction of the values of the couplings, however, 
becomes more involved than in leading order. 

3 Diagonalisation in the couplings 

We shall now propose a method to analyse the data which presents several 
advantages in view of the basic problem posed by the large number of un- 
known three-boson couplings: with limited event statistics significant error 
bounds can only be obtained for subsets of the coupling parameters, but im- 
posing constraints on the couplings to reduce their number entails a loss of 
information that cannot be retrieved. In view of this it should be advanta- 
geous to use a parametrisation of the couplings which in a given process and 
at a given cm. energy has the following properties: 

1. It allows to find observables which are only sensitive to one particular 
coupling parameter. 

2. The induced errors on the couplings determined from these observables 
are statistically independent. 

With this we can on the one hand give single errors for each parameter, on the 
other hand we can recover from the single errors the multidimensional error 
of the full set of couplings, having avoided the loss of information incurred 
by imposing constraints. From the single errors we can also directly see 
which combinations of couplings in more conventional parametrisations can 
be measured with good accuracy and to which one is rather insensitive. 

Let us remark that in the leading order analysis there is a set of observ- 
ables satisfying point 1. in any parametrisation of the couplings. The linear 
combinations 

C, = Y.c-\,0, (17) 



of our optimal observables (13) are only sensitive to gi for each i (cf. also 



P^). The errors on the couplings determined from these observables are, 
however, in general not uncorrelated; in fact their correlations are the same 
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(a) 



(b) 






Figure 1: Example of the l-o" ellipsis for two extracted parameters in a 
parametrisation where their errors are correlated (a) or uncorrelated (6). 
The single errors on the couplings are given by the projection of the ellipses 
on the coordinate axes. 



as those obtained with the original set Ot. This can be seen as follows: going 
from the Oi to the Cj we must replace 



6. 



t] ! 



V{0) -^ V{C) = c-^-V{0)-{c-^) 



l\T 



SO that we have from fO 



Vigy 



o 



Nc' ■V{0)-^-c = NV{C) 



-1 



V{9) 



-1 



(19) 



In such a case the single errors give an incomplete picture of the situation if 
correlations are large. This is illustrated in fig. |I] (a), where the 1-cr ellipsis 
for two parameters is shown. Their single errors are given by its projection 
on the coordinate axes and in our example are both rather large. Some linear 
combinations of them are however measurable with much better precision, 
which one can only recognise if both errors and their correlations are given. 
In fig. [1| (6) where a set of couplings leading to uncorrelated errors is used 
the situation is much simpler. Note also that the number of correlations, i.e. 
off-diagonals in V{g), is yet modest for two couplings but increases rapidly 
with their number. 

We will now first see that a parametrisation of the couplings satisfying 
both points 1. and 2. above can be found in idealised circumstances, and then 
mention the restrictions one will encounter under more realistic assumptions. 
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If the leading order analysis is a sufficiently good approximation the so- 
lution to our problem is easily found. Starting from a set of couplings gi and 
the corresponding optimal observables Oi in ( ]T3| ) we can go to another set g[ 

by 

g' = A-ig , (20) 

where we use vector and matrix notation. The coefficients in the expansion 
of the differential cross section and the optimal observables transform as 
follows: 

s; = A^Si 

O' = A^O (21) 

Let now Oi be an arbitrary set of observables related to gi and define the 
corresponding 0[ related to g[ as in (0). Then we have for the matrices 
relevant for our analysis the following transformation properties: 



c' 


= A^-c-A 


V{0') 


= J^-V{0)-A 


vr' 


= A^ . V{g)-' ■ A . 



(22) 

As shown in [^ our optimal observables satisfy V{0) = c and V{g)^^ = Nc 
so that for them one can choose a transformation A which diagonalises all 
three matrices. This new set g' of parameters obviously has the properties 
1. and 2. we were looking for. 

Beyond the linear approximation of @ the expectation value of O'^ will 
still receive contributions from several couplings. In fact there is no set of 
observables for which the full nonlinear expression in (^) satisfies point 1. 
exactly, because the denominator involves quadratic terms in all couplings, 
and this cannot be changed by any linear transformation of the couplings. If 
on the other hand the statistical errors are too large the covariance matrix 
V{g) will not give a good picture of the errors as we discussed in sec. ||, and its 
diagonalisation will not ensure point 2. In the case however where nonlinear 
effects in the determination of the couplings and their errors are not too large, 
i.e. where the leading order expressions are a good ffist approximation both 
points 1. and 2. above will still be approximately satisfied in a full nonlinear 
analysis. 

We remark that if one has some previous estimates gi of the couplings 
that considerably deviate from zero one may reduce nonlinear effects in the 
determination of the gi by working with an expansion of da/dcj) around the 
cji as shown in sec. |^; in our diagonalisation programme one will then use 
couplings gi — gi instead of gi, the matrix c instead of c etc. 
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To the extent that the observables (p!3| ) are constructed from expressions 
of 5*0 and Sij which are only approximations of those that determine the 
experimentally observed kinematical distribution the matrices c, V{0) and 
N~^V{g)~^ will not quite be the same and cannot be diagonalised at the 
same time. One can then diagonalise either V{0) or V{g)~^ because they 
are by definition symmetric and positive definite, whereas c is not necessarily 
so. Again, unless such effects are large one will end up with a matrix c' that 
is not diagonal but has relatively small off-diagonals. 

It should also be borne in mind that the covariance matrix V{g) only gives 
the statistical errors on the couplings, so that even if it is exactly diagonal 
the final errors may be correlated due to systematics. 

3.1 Simultaneous diagonalisation of the correlation 

matrix and the quadratic term in the total cross 
section 

The choice of transformation in (^) to (^) is not unique if one does not 
require A to be orthogonal.^ We see in fact no strong argument in favour 
of an orthogonal transformation and remark that the various parametrisa- 
tions of the WW'j and WWZ couplings in the literature are related by non- 
orthogonal linear transformations. The freedom to choose A can be used to 
impose additional conditions on the transformation, and the one we propose 
here is that the transformed quadratic coefficient 

a'^ = A^ ■a2-A (23) 

in the integrated cross section be the unit matrix. In terms of the new 
couplings one then has 

8 28 

a/ao = l + E^M^: + E(^:)' , (24) 

where we choose the numbering such that g'l to g'g belong to symmetry class 
(a) introduced in sec. p?T| i.e. they are the CP and CPT even couplings. 
Only these appear linearly in the cross section, whereas all couplings give 
a quadratic contribution with coefficient one. Having (J2 j,- = Sij leads to a 
convenient simplification of (^), (0). Moreover, the measurement of the total 
cross section gives complementary information on the unknown couplings. 



(25) 



Rewriting 


(11) as 


















a/ao = 


= 1- 


8 

->; 


4 


8 

+ >; 


;.;* 


2 


2 28 

+ > W' 
i=9 



^Orthogonal transformations have been used in the second paper of ref. M. 

13 



we see that measuring a cross section a^xp within an error Aa constrains the 
couphngs to be in a shell between two hyperspheres with centre at 

g[ = -^ (2 = 1,. ..,8), g[ = Q (z = 9,...,28) (26) 

in the space of all couplings as shown in fig. 0. Their radii are given by 

r> = — , r< = — ^ . (27) 

O"0 CTq 



Here 



^min '-^Q 



1- 



i^'u 



8 (;:,' \2^ 

E 

1=1 



:28) 



is the smallest value the cross section can attain; that such a minimum exists 
has been pointed out in |]T5|. If in (^) r| is positive but r\ negative the 
couplings are inside the hypersphere with radius r>, and if both r> and r| 
are negative the ansatz ([l|) for the cross section is inconsistent with the data 
within the error Act. 

Such constraints can be useful to find the physical set of couplings when 
the solution of (^ from the measurement of the optimal observables is not 
unique. If they are strong enough they might even restrict the couplings to 
the region where (^ can be linearised and thus simplify their extraction. One 
can of course use the information from the integrated cross section working 
with any set of couplings, but again the situation is particularly simple with 



the form (p^). 

We note that the information from the total rate is complementary to 
what is extracted from the mean values of our observables, which involve 
normalised kinematical distributions. From an experimental point of view 
their respective measurements will presumably have quite different system- 
atic errors. Let us also recall that the the measurement of the mean values 
Oi times the number A^ of events obtained with a fixed integrated luminosity 
combines the information of both [Q . A nonlinear data analysis as presented 
in sec. ^ can also be done in this case. We draw attention to the fact that 
unphysical solutions of equation (|^) for Oi and of its analogue for NOi will 
in general be different. We shall however not elaborate on this point here. 

Another aspect of the couplings with the property (|^) is the following. It 
is well known that constant coupling parameters deviating from the standard 
model tree level values lead to amplitudes that violate unitarity |T^. The 
coefficients ai^j and (T2,jj in the total cross section a increase strongly with 
the e"'"e~ cm. energy a/s and the couplings Qi must vanish as s becomes large 
to ensure a decent high-energy behaviour of a. In our new parametrisation 
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Figure 2: The measurement of the integrated cross section constrains the 
couphngs to a shell between two hyperspheres with radii r> and r< given by 



(|27|). g'i belongs to symmetry class (a) introduced in sec. ^]T] and g'j to class 
(b), (c) or (d). In the example shown here the measurement is compatible 
with both couplings being zero. 



the quadratic coefficients in (jj are energy independent, and in this sense the 
new couplings are at a "natural scale" at every energy. 

To complete this section we show that a transformation with the proper- 
ties we require always exists, i.e. that we can find a matrix A that diagonalises 
V{g)~^ in ( P^D and transforms (J2 in ( PB| ) to the unit matrix. The argument 
is analogous if one replaces V{g)^^ with V{0). By construction both V{g)~^ 
and o"2 are symmetric and positive definite, so our problem is the same as 
finding normal coordinates for a multidimensional harmonic oscillator in clas- 
sical mechanics (cf. e.g. |[ll|]). To make this analogy transparent let us write 



(29) 



T = (5"2 and V = V{g) , we then have to find A so that 



A^-T-A 



A^ -V- A = D 



with D being diagonal. The elements di of D are generalised eigenvalues of 
V satisfying 



V a,; = di TsLi 



(30) 



15 



where a^ is the i-th column vector of A. The solution is well known to be 

A = T-i/2 . u (31) 

where U is the orthogonal matrix that transforms T~^/^ ■ V ■ T^^/^ to D. Of 
course one need not use ( |3TD in practice as there are convenient algorithms 
available to find A and D. In our numerical calculations we have used the 
routine Eigenvals of the algebraic package MAPLE. 

4 Numerical examples 

We will now give some numerical examples of our method described in the 
previous section. In this section we will stay within the framework of a 
leading order analysis of the observables. We start from the results in |^, 
where the sensitivity of optimal observables for semileptonic ly py-decays was 
calculated. We assume a full kinematical reconstruction of the final state, 
except for the ambiguity one is left with if the jet charge is not known. For 
the standard model cross section we use the Born approximation and neglect 
effects of the finite W width. To describe the triple boson couplings we take 
the form factors /7, ff [i = 1 . . .7) oi ^] deviations from their standard 
model tree level values will be referred to as "anomalous couplings" . 

Let us first look at a cm. energy of ^/s = 190 GeV, which will be at- 
tained at LEP2. The coefficient matrix c can be found in table 4 of |^ and 
in table ^ here we give the diagonal elements of the transformed matrix c', 
ordered according to the symmetry of the corresponding observables. The 
one standard deviation ellipsoid is diagonal in the couplings g'^, thus its in- 
tersections with the coordinate axes equal its projections on these axes. The 
errors dg'i setting all other g', to zero are then equal to the errors Ag'^^ where 
all other g'j are arbitrary. They are given by 

^9[ = , ^ = -^ (32) 

^V{g')-\i V^ 

and are listed in table |^ for an integrated luminosity of 500 pb^^. 

We immediately remark that a negative diagonal element occurs in the 
transformed coefficient matrix, which is not allowed because c' = V{0') is 
a covariance matrix and thus positive definite. We encounter here a prob- 
lem of numerical instability: small errors in the calculation of the original 
matrices c and o"2 can have a large effect on the smallest generalised eigen- 
values c^j and their eigenvectors, even to the point that eigenvalues come 
out with the wrong sign. This is not only a problem of our particular way 

16 



of diagonalisation, but also occurs if one diagonalises c with an orthogonal 
matrix; we find that one of the usual eigenvalues of c in the subspace of 
couplings with symmetry (6) is negative. Such instabilities can cause large 
errors in the matrix inversion of c and V{0). One needs V{0)~^ to calculate 
the error on the extracted couplings as can be seen from @ and (0), and 
large errors on c~^ can lead to large uncertainties in the extracted couplings, 
irrespective of whether c~^ is explicitly used to solve the system (|^. One 
will of course aim to calculate c and V{0) with best possible precision, but 
such an effort has limits, in particular if they are determined from simulated 
events and include for instance radiative corrections or detector effects. On 
a more fundamental level any calculation of these matrices will only be an 
approximation of the "exact" ones that correspond to the kinematical dis- 
tributions seen in experiment. In this sense it seems quite inevitable that 
small eigenvalues (the usual or our generalised ones) of c and V{0) and their 
eigenvectors are sensitive to imprecisions in the calculation and can lead to 
large errors or uncertainties in the data analysis. This holds of course even if 
one does not obtain eigenvalues with the wrong sign. We think that also in 
view of this a diagonalisation is useful, not because it solves the problem but 
because it makes it explicit! It allows to easily identify those combinations 
of couplings which have small corresponding eigenvalues in c and V{0) and 
will be the most unsafe ones in the analysis. From (^) we see that they are 
those combinations for which the statistical errors will be largest. Here the 
most unsafe coupling parameter is g'^Q. One might thus choose to exclude it, 
and possibly other couplings, from the analysis and work in the remaining 
subspace of the g'^ where the numerics is more stable and where in any case 
the experiment is most sensitive. We will come back to this in sec. ^. 

In tables |l| and |^ we find that the range of sensitivities is quite large, 
typically spanning several orders of magnitude in the c^j. We can actually 
identify the role of those form factors whose c^j are small. To this end we 
pass from the usual form factors //, ff to the combinations which appear 
in the amplitudes for e~^e~ -^ W^W~ with left or right handed electron 
polarisation, respectively: 

//- = 4sin2^H'/7 + (2-4sin2Me/f, 

/f = 4sin2^H^/7-4sin2^H/e/f, (33) 

where ^ = s/{s — M|) and 6^ is the weak mixing angle. The matrix o"2 
is block diagonal in the left and right handed form factors because different 
electron polarisations do not interfere in the cross section, but the coefficient 
matrix c is not. If one sets all right (left) handed anomalous couplings to zero 
then one has to diagonalise c in the left (right) handed subspace. Doing this 
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Table 1: Diagonal elements c^^ of the transformed coefficient matrix at ^/s = 
190 GeV. They are ordered by the symmetry classes (a) to [d) introduced in 
sec. pri| ; the first row contains c[ i to Cgg, the second Cgg to c'^g ig etc. Note 



that the c^j are by construction positive; the negative value for c'^g ^g is due to 
numerical integration errors in the original matrix c. This problem is further 
discussed in the text. A similar comment applies to tables |^, ^ and ^. 



(a) 


1.5 


1.4 


0.83 


0.70 


0.32 


0.10 


0.027 


0.019 


(b) 


1.1 


1.0 


0.81 


0.68 


0.028 


0.0093 


0.0056 


-0.00045 


(c) 


0.80 


0.67 


0.28 


0.028 


0.013 


0.012 






id) 


1.4 


1.2 


0.85 


0.13 


0.039 


0.017 







Table 2: 1-cr errors Ag^ on the extracted couplings corresponding to the 
coefficients c^j in table |1|. They are calculated from (^) for an integrated 
luminosity of 500 pb~^. 



(a) 


0.015 


0.016 


0.020 


0.022 


0.032 


0.058 


0.11 0.13 


(b) 


0.017 


0.018 


0.020 


0.022 


0.11 


0.19 


0.24 — 


(c) 


0.020 


0.022 


0.035 


0.11 


0.16 


0.17 




id) 


0.015 


0.017 


0.020 


0.051 


0.093 


0.14 





Table 3: Diagonal elements c^j of the coefficient matrix restricted to the left 
or right handed subspace of the couphngs as explained in the text. The values 
in the left handed subspace differ from the corresponding ones in table 1 by 
at most 3%. 







left handed 






right handed 




(a) 


1.5 


1.4 0.83 


0.69 


0.35 


0.12 0.030 


0.019 


(b) 


1.1 


1.0 0.81 


0.68 


0.030 


0.23 0.11 - 


-0.00042 


(c) 


0.80 


0.67 0.28 




0.030 


0.018 0.012 




id) 


1.4 


1.2 0.85 




0.13 


0.042 0.018 
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we obtain the results shown in table 0. In each of the classes (a) to {d) we 
find a clear correspondence of the generalised eigenvalues in the left (right) 
handed subspace with the largest (smallest) ones of the full matrix c given in 
table [^. The form factors g'^ to which one is most sensitive thus correspond 
predominantly to left handed combinations, whereas the right handed com- 
binations are more difficult to measure. This confirms our findings in and 
can be explained by the missing neutrino exchange graph for right handed 
electrons which can give a large interference with anomalous triple boson cou- 
plings. As a word of caution we remark that the values given in table ^ do 
not correspond to those for left or right handed electron beams, because for 
unpolarised beams both electron helicities contribute to the standard model 
cross section even if anomalous couplings corresponding to one helicity are 
(assumed to be) zero. We found however in p that the difference between 
c for a left handed electron beam and the left handed submatrix of c for 
unpolarised beams is small, again because the right handed standard model 
contribution is small compared to the left handed one due to the missing 
neutrino exchange graph. 

Let us now come to the integrated cross section. We first give it in the 
parametrisation by // and ff, where the standard model tree values are 
fi=fi=^if2 — fi — '^ ^^d zero for all other form factors: 

cr/do = 1 

+ 0.022 (Re/7 - 1) + 0.013 (Re/f - 1) - 0.031 Re/2'' - 0.010 Re/f 
- 0.074 (Re/3^ - 2) - 0.019 (Re/3^ - 2) + 0.0058 Re/5'' + 0.0092 Re/5^ 
+ {quadratic terms} (34) 

We do not give the full matrix (J2 for the quadratic terms here, but remark 
that the absolute values of its elements are between 0.001 and 0.3 and its 
eigenvalues between 0.004 and 0.4. In our new parametrisation we have 

a/ao = 1 + 0.18^1 + 0.16 ^2 -0.0053 ^3 -0.052 c/; 

28 

- 0.15 g'r^ + 0.071 g'^ - 0.029 g', + 0.0091 g'^ + Y.{g[f 

i=l 

= 1-0.023 

+ {g[ + 0.092)2 ^ (^^^ ^ 0.082)2 ^ ^^^^ _ 0.0026)^ + (g'^ - 0.026)2 
+ {g', - 0.073)2 ^ ^^, ^ Q_Q35)2 ^ ^^, _ q_q^4)2 ^ (^/ ^ 0.0046)2 

28 

+ E(^0'- (35) 

Comparing with the first row of table |l| we see that couplings whose lin- 
ear contribution to the cross section is relatively small can give a relatively 
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large linear contribution to their optimal observable and vice versa. A mea- 
surement with an integrated luminosity of 500 pb^ at y/s = 190 GeV that 
finds the cross section equal to its standard model (Born level) value ctq will 
have a relative statistical error Aa/aexp = 1/VN = 0.018. According to 



(|27|) this measurement would constrain the couplings g^ to be between hy- 
perspheres with radii r< = 0.066 and r> = 0.202 in their 2 8- dimensional 
space. Comparing with table ^ we see that the thickness r^ — r^ = 0.136 of 
this shell is of the same order of magnitude as the largest statistical errors 
to be achieved with our optimal observables, which give the extensions of 
the 28-dimensional error ellipsoid in a linearised analysis. This reflects the 
well-known fact that the integrated cross section is clearly not as sensitive 
to anomalous couplings as the detailed kinematical distributions whose in- 
formation we extract with our observables. It gives however an additional 
constraint that should provide a means to constrain the couplings such as 
g[Q which are practically unmeasurable by normalised observables. Also we 
should get a useful cross check and help selecting the physical solution of (|^) 
in a full nonlinear analysis. 

Let us now apply our method to a typical LC energy y/s = 500 GeV with 
an integrated luminosity of 10 fb~^ (cf. tables ^ to |^). Here we encounter 
the particularity that the coefficients of the form factors f2' and /^ ' in 
the scattering amplitude grow faster with energy than those of the other 
couplings [l|] . As a result the off-diagonal matrix elements in c between /g^' 
and another coupling are about a factor of 10 larger than matrix elements 
not involving f2' , and elements involving only /^' are larger by yet another 
factor of 10. The (usual) eigenvalues of c in the CP conserving sector span 6 
orders of magnitude. In the CP violating sector the situation is less dramatic, 
but still elements of c involving /y' alone are about a factor of 10 larger than 
the others. The same phenomenon is found in the cross section, which reads 

a/ao = 1 

+ 0.45 (Re/7 - 1) + 0.23 (Re/f - 1) - 8.3 Re/2^ - 4.3 Re/f 

- 0.58 (Re/3^ - 2) - 0.31 (Re/3^ - 2) + 0.056 Re/5'^ + 0.070 Re/5^ 

+ {quadratic terms} (36) 

where (T2 has elements with absolute values between 0.04 and 1700 and eigen- 
values between 0.02 and 1800. After our simultaneous diagonalisation the 
range of the matrix elements c[^ is significantly smaller as can be seen in 
table ^. The cross section now reads 

a/ao = 1-0.31^; + 0.057^2 + 0.0018^3 -0.014c/; 

- 0.026 (7^ - OMAg', - O.OUg', - 0.013 .7^ + J2{g: 



28 

')' 
1=1 
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= 1-0.026 

+ {g[ - 0.16)2 + {g'^ + 0.028)2 + (g'^ + 0.00091)^ + (g'^ - 0.0068)^ 
+ (^5 - 0.013)2 + (g'Q - 0.022)2 ^ ^^^/^ _ 0.0069)2 + (g'^ - 0.0066)2 

28 

+ T.i9'f . (37) 

1=9 

Of course, the g'^ are now in general energy- dependent for constant form 
factors /7' . As we mentioned in sec. p.l| the rapid growth with s of the 
coefficients of anomalous parts of the couplings //' in the amplitude has 
been "absorbed" into the coupling parameters by the transformation of o"2 to 
the unit matrix. As a result the linear coefficients in the cross section and the 
elements of the coefficient matrix c at y/s = 500 GeV have the same order of 
magnitude as at LEP2 energies, and the large differences in the s-dependence 
of the coefficients for different couplings have been evened out. 

A measurement at the LC that gives the standard model cross section 
constrains the couplings to be between 28-dimensional hyperspheres with 
radii r< = 0.140 and r> = 0.182, taking the relative statistical error of 
0.0067 on the cross section one would obtain with an integrated luminosity 
of 10 fb^ . The thickness of the shell, r^ — r^ = 0.042, is again of the 
order of the largest statistical errors on the g^ one can obtain with optimal 
observables (cf. table ^. Moreover such a small statistical error is likely 
to be small compared with systematic errors, so that the sensitivity of the 
integrated cross section will even be less. 

Finally we remark that like in the case for ^/s = 190 GeV those couplings 
g'j^ which give the largest statistical errors in the optimal observable analysis 
are predominantly related to right handed combinations of form factors as 
can be seen from the comparison of tables ^ and ^. 

5 Simultaneous diagonalisation in practice 

Let us sketch how our method of simultaneous diagonalisation might be used 
in practice. 

1. One first has to choose which matrix to diagonalise simultaneously with 
(T2. These matrices need not be the same ones to be used in the data 
analysis itself but may be calculated under further approximations. 
Covariance matrices for the observables and extracted couplings can be 
evaluated for zero gi as our entire procedure will only have its desired 
properties if nonlinear effects are not too large. If one uses the same 
approximation of the differential cross section (|^) for the construction 
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Table 4: As table |l|, but for ^/^ = 500 GeV. 



(a) 


1.4 


1.2 


0.74 


0.65 


0.33 


0.11 


0.056 


0.033 


(b) 


1.3 


1.0 


0.79 


0.29 


0.097 


0.056 


0.0092 


-0.0013 


(c) 


1.2 


0.58 


0.32 


0.066 


0.027 


0.013 






id) 


1.4 


1.0 


0.83 


0.22 


0.033 


0.025 







Table 5: 1-cr errors Ag'^ on the extracted couplings, corresponding to the 

III I.M I HI-' lOJ 1(11 Mil I I I l.f-'U IM l.f-'l 1 I I I I I I I I 11 IKI lA/ 1 II III 11^ 



coefficients c[^ in table H for an integrated luminosity of 10 fb 



(a) 


0.0056 


0.0062 


0.0078 


0.0083 


0.012 


0.020 


0.029 0.037 


(b) 


0.0060 


0.0066 


0.0075 


0.013 


0.022 


0.028 


0.070 — 


(c) 


0.0062 


0.0088 


0.012 


0.026 


0.041 


0.059 




id) 


0.0057 


0.0067 


0.0074 


0.014 


0.037 


0.043 





Table 6: As table |^ but for y/s = 500 GeV, to be compared with table 







left handed 




right handed 




(a) 


1.4 


1.1 0.72 0.63 


0.50 


0.17 0.062 


0.044 


ib) 


1.3 


1.0 0.79 0.26 


0.11 


0.083 0.023 


-0.0012 


(c) 


1.2 


0.58 0.32 


0.076 


0.031 0.013 




id) 


1.4 


1.0 0.83 


0.24 


0.040 0.026 
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of the optimal observables (|T3|) and the calculation of c, V{0) and 
N^^V{g)~^ then the latter are all equal and can be diagonalised at the 
same time. Otherwise one has to choose a positive definite symmetric 
matrix for the diagonalisation, i.e. one of the covariance matrices. The 
calculation of V{g) or of its inverse from (12) involves however a matrix 



inversion and might suffer from numerical instabilities, so presumably 
the best choice will be V{0). 

2. In the next step one carries out the simultaneous diagonalisation of 
the chosen matrix and (J2 as described in sec. p.l| and determines the 
transformation matrix A in ( pO] ) to p^). At this point it will be useful 
to test the numerical stability of the transformed matrices, for instance 
by re-calculating them in the new basis of couplings or by repeating 
the diagonalisation procedure with slightly modified initial matrices. 
One might choose to discard some of the new couplings g'^^ and the 
corresponding observables from the analysis if the corresponding matrix 
elements are found to be instable. This does not mean that one has 
to set these couplings to zero. From the measurement of the total 
cross section one will obtain limits on them, which will also allow to 
control the contribution they can give to the mean values of those 
observables that are kept in the analysis because the matrix c' is not 
exactly diagonal and because of nonlinear terms in (|^. 

3. In the new parametrisation of the couplings one then carries out the 
analysis of the data. Here V{0'), c', o"2 and the other coefficients in 
(^ will be determined under the most realistic assumptions and with 
the best precision one can afford. They will not be exactly diagonal in 
practice, but should have small off-diagonal elements if the approxima- 
tions made in step 1. and in the construction of the optimal observables 
are sufficiently good. 

4. One can then give both single and multidimensional errors on the mea- 
sured coupling parameters g^. At this stage one can also present the 
results in other, more conventional parametrisations of the couplings 
and in particular compare with the measurements of other experiments, 
restricting oneself to whatever subspace of couplings might have been 
chosen there. 
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6 Summary 

In the first part of tliis paper we liave sliown liow to extract coupling param- 
eters from tlie measured mean values Oi of appropriate observables without 
the approximation that the couplings are small. Errors on the couplings can 
be obtained from a x^"fit of the Oi. If one puts constraints on the cou- 
plings in order to reduce their number the method also gives an indication of 
how compatible these constraints are with the data. The "optimal observ- 
ables" discussed in have statistical errors equal to the smallest possible 
ones to leading order in the coupling parameters gi. Beyond the leading 
order approximation one can obtain more sensitive observables if one has 
some previous estimate Qi for the couplings, expanding the differential cross 
section around g^ instead of zero and constructing observables from the cor- 
responding expansion coefficients. In the appendix we show that up to linear 
reparametrisations the choice of optimal observables is unique: any other set 
of observables must give bigger (statistical) errors. 

In a second part we have proposed to perform the data analysis using a 
particular parametrisation g^ of the couplings, which is specific to the process 
and its cm. energy. It is obtained from the initial set g^ by a linear trans- 
formation which diagonalises the covariance matrix V{0) of the observables 
and transforms the matrix o"2 of quadratic coefficients in the integrated cross 
section (Q) to unity. In an idealised framework each optimal observable O'^ 
for this parametrisation is only sensitive to one coupling, and the statistical 
errors on the extracted couplings are uncorrelated. Under realistic circum- 
stances both properties can be approximately satisfied provided that the 
analysis stays in a region of parameter space where the dependence of the 
mean values O'^ on the couplings is not far from linear. Various matrices are 
then approximately diagonal which should generally facilitate the data anal- 
ysis. In particular one can directly give errors on single or a small number 
of couplings, which will be necessary to obtain statistically significant results 
with a limited number of events. At the same time one can readily present 
multidimensional errors in parameter space, which is essential to compare 
with the results of measurements that impose various different constraints 
on the couplings. Having approximately diagonal matrices also allows to eas- 
ily identify those directions in parameter space which can be measured best 
and those for which the statistical errors will be large and which are likely to 
be associated with numerical instabilities, for example in matrix inversions. 
One can thus recognise and seek to remedy such problems in an early stage 
of the analysis. 

The measurement of the total cross section a gives valuable complemen- 
tary information on the coupling parameters. Its dependence on the cou- 
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plings is particularly simple in the parametrisation we propose since the 
quadratic contributions are {g'^)^ times the standard model cross section ctq, 
i.e. they have the same form for all couplings. A measurement of a will then 
restrict the g[ to a shell between two hyperspheres in parameter space. 

We have given some numerical examples of our method applied to the 
semileptonic decay channels in e~^e~ -^ W^W~ . In particular we find that 
the couplings g[ which can be measured best with unpolarised beams predom- 
inantly appear in the amplitude for left handed electrons (or right handed 
positrons), and that the g[ with the largest statistical errors mainly corre- 
spond to the opposite lepton helicity Comparing our results at LEP2 and 
LC energies we see that the coefficients in the linear contributions of the 
couplings g[ to our observables and to the integrated cross section change 
much less with energy than in usual parametrisations. This is because in 
the new parametrisation the quadratic coefficients in the normalised cross 
section a /gq are by construction energy independent. 
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Appendix: Uniqueness of optimal observables 



In this appendix we show that the set of observables ([Iq) , obtained from 
expanding the differential cross section about the actual values of the cou- 
plings, is unique in the sense that up to the hnear reparametrisations (^) it 
is the only set of n integrated observables which in the limit of large N leads 
to the minimum error on the n extracted parameters. 

To keep our notation simple we give the proof for the case that the actual 



25 



values of the Qi are zero. The expectation value and covariance of functions 
/(0) and g{(f)) are then given by 



In the general case one has instead of So the zeroth order coefficient "So (14) 
from the expansion about the appropriate values Qi. 

For large N the covariance matrix for the extracted couplings is given by 
(|T^). Under a linear reparametrisation of observables, 

O,(0) ^ OK0) = E (^^,^A<P) + h ' (39) 

i 

where a^j and 6j are constants and the matrix a^j is nonsingular, the matrices 
c from (^ and V^(0) transform according to 

c' = a ■ c ^ 
V{0') = a ■ V{0) ■ J . (40) 

From (|l^) we see that the covariance matrix V(^g) is unchanged under such a 
transformation. For our proof we can hence restrict ourselves to observables 
with mean value 

Eo[0,] = (41) 

and with a coefficient matrix Cjj = bij. From (^ we then have the condition 

K)[0., ^i,,/^o] = <5., (42) 

and the error on the extracted couplings is given by 

V{g),, = N-'V{0),, = iV- Vo[0„ O,] . (43) 

From [^ we know that the optimal observables (0) lead to the smallest 
possible error on the Qi, given by the Cramer- Rao bound. To satisfy our 
conditions (41) and (^21) we take the linear combinations 



M<P) = T.i^optr\j (5i,,(0)/5o(0) - Eo[S^,/So]) (44) 

j 

with 

{copt)ij = Vo[Si^i/So, Sij/So] . (45) 

We assume that the functions Si^i are linearly independent. Otherwise some 
of the parameters Qi are superfluous and can be eliminated; our assumption 
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is thus that the gi are an independent set of parameters for the anomalous 
couphngs. Linear independence of the 5'i^j(0) guarantees that Copt is nonsin- 
gular, which has tacitly been used at several instances in our paper. The set 
Vi is related to the optimal observables Si^i/ Sq by a linear transformation 
(|39|) and thus gives the same optimal error matrix V{g). 

The covariance Vq [/,(?] defines a scalar product on the Hilbert space Ti 
of sufficiently smooth functions of with the property -Eo[/] = 0.| The 
functions Vi span a subspace Ti^ of Ti, and we define Ti^^ as the orthogonal 
complement of 7i^ with respect to the scalar product Vq [/,(?]. Any set of n 
observables satisfying (O) can then be written as 



O, = Ol + Of (46) 

with ol G 7i^, Of G Ti^^ . Further decomposing (9/ = J^jO'ijI^j and using 
the constraint (42) we obtain aij = 6ij, i.e. 

O- = V, . (47) 



Finally, we have from (|43|) 



V{g),, = N~Wo[V„V,]+N-Wo[0'/,Of] . (48) 

The first term gives the error on the couplings for the optimal observables Vi , 
which is minimal. If the observables Oi have minimal error, too, the second 
term must be zero, so that for each i we have Vo[Cf , Cf ] = and thus 

Of = , (49) 

which completes our proof. 

In sec. ^TT| we mentioned that instead of Oi one may use the product 
NOi measured with fixed luminosity to extract the couplings 0. By an 
argument analogous to the one of this appendix one finds that up to linear 
reparametrisations our observables ([T6|) are again the only optimal ones. In 
this case linear reparametrisations have to be homogeneous, i.e. one must 
have 6j = in (0), since adding constants to the observables can change the 
induced errors on the coupling parameters. 
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